Bilayer Rydberg atoms as a quantum simulator for unconventional superconductors. 
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We show how systems of cold fermionic Rydberg atoms in a bilayer lattice can be used as a 
quantum simulator for electron-phonon interactions in the presence of strong electronic correlation: 
a scenario found in many unconventional superconductors. We discuss the experimental system 
for making such a simulator, show that the resulting Hamiltonian can be mapped to an extended 
Hubbard-Holstein model, and make comparison with numerics. Fermions can be simulated, and 
the proposed quantum simulator can cope with systems where numerics are difficult to apply. 



Besides the cuprates, there are several superconductors 
with high transition temperatures, many of which have 
important roles for electron-phonon interactions, and re- 
pulsion driven correlated electron phenomena such as an- 
tiferromagnetism. Fulleride superconductors of the fam- 
ily A3C60 have phonon driven transition temperatures 
PQ of up to 40K [2], but also exhibit antiferromagnetism 
at appropriate dopings and structures [3]. Bismuthate 
superconductors with T c > 30K [JJ have electrons well 
localized to atomic orbitals with a low density of states 
at the Fermi surface, with superconductivity probably 
due to strong couplings to the lattice [5j (as evidenced 
by a large isotope shift High transition tempera- 

tures have also been achieved in the borocarbides [5J |5] 
(T c = 23K) and chloronitrides [TU] (T c = 25K) which 
are also thought to be driven by electron-phonon cou- 
pling. Even the cuprate superconductors, where super- 
conductivity is thought by many to be driven by anti- 
ferromagnetic fluctuations [IT] , show isotope shifts [H] 
and other effects such as kinks [13] that may be at- 
tributed to some kind of non-trivial interplay between 
strong correlations and lattice vibrations. 

Owing to the importance of electron-phonon interac- 
tions in condensed matter, reliable numerical methods 
have been sought, but even very simplified models |14j 
are extremely hard to simulate. The difficulty lies in 
dealing with the potentially infinite number of phonons 
associated with even a single electron. Another key prob- 
lem in many advanced materials is that electrons are 
localized to atomic orbitals and accordingly the Fermi 
energy is small. This localization means that dimension- 
less electron-phonon couplings are relatively strong and 
phonon frequencies can be large at around 10% of the 
Fermi energy. Even moderate couplings with interme- 
diate frequency phonons can lead to consequences that 
cannot be predicted with perturbation theory and other 
analytics. Such couplings can cause additional difficul- 
ties for numerics that often work most efficiently in the 
extreme limits of weak or strong coupling. 

Cold atom quantum simulators offer an important new 
approach to the study of correlated electron phenomena 
without the limitations of computational or analytical 
techniques. For example, cold atoms have recently been 



used as quantum simulators to investigate standard mod- 
els of condensed matter physics such as the Hubbard 
model of strong local Coulomb repulsion [T^l D2>] ■ The 
use of cold atoms to simulate Hubbard models has led to 
direct observation of important phenomena such as the 
superfluid to Mott insulator transition [T71 [TH] . 

Quantum simulation of lattice effects has proved dif- 
ficult to implement. The existence of quantum simu- 
lators with a high degree of control over the form of 
interactions has the potential to provide significant in- 
sight into the subtle interplay between electronic correla- 
tion, lattice vibration and phenomena such as supercon- 
ductivity and colossal magnetoresistance. In this letter, 
we propose an approach to simulating fermionic Hub- 
bard models extended to include strong long-range in- 
teractions with lattice vibrations, by discussing Rydberg 
states of cold atoms in bilayer lattices. Such a simu- 
lator is likely to shed light on a wide range of uncon- 
ventional superconductors with transition temperatures 
greater than 30K, and might resolve some aspects of on- 
going debates on cuprate mechanisms. The simulator 
is introduced and the mapping to a corresponding ex- 
tended fermionic Hubbard-Holstein model is carried out. 
We perform simulations to show how the physics of the 
system is comparable with models of high temperature 
superconductors. We then discuss how to read proper- 
ties of the resulting correlated fermion problems from the 
simulator. 

Several schemes for quantum simulation of interactions 
between electrons and phonons have been suggested. A 
proposal to bathe an optical lattice in a BEC [HI [20] 
has led to the observation of polaron effects [31]. In- 
teractions with high energy phonon states of Rydberg 
ions are used as part of the mapping in proposed simu- 
lators for spin systems [22, 23J. Li and Lesanovsky have 
discussed structural distortions associated with exciting 
high energy Rydberg states in cold ion crystals [24] . Ry- 
dberg atoms have been proposed as a way of simulating 
polaron effects in strongly deformable materials [35] . The 
use of cold polar molecules to obtain Holstein polaron ef- 
fects has been discussed [26 . A scheme for simulation of 
the Frcnkcl- Kontorova model to probe friction and en- 
ergy transport [37] has also been proposed. The scheme 
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FIG. 1. System of bilayer Rydberg cold atoms for simulation 
of strong correlations and interactions between fermions and 
phonons, annotated with Hamiltonian terms. 



proposed in this letter goes beyond these as it is capa- 
ble of simulating interacting Fermions, rather than the 
interacting Bosons in the earlier schemes. 

We begin by discussing how electron-phonon interac- 
tions can be simulated in a system of cold, highly excited 
(Rydberg) atoms in a bilayer lattice. We assume that 
the Rydberg atoms interact only via dipolc-dipolc inter- 
actions of the form: 

Vki = - In /i2 „„ , (1) 



|i4--R;| 3 ' 

where /x is the dipole moment on the Rydberg atom and 
Rk is the vector to the fcth atom in the lattice. The dipole 
moment may be written in terms of the coefficient C3 as 
/i = v2C 3 , and depends upon the Rydberg states chosen 
for the experiment. The ground state atoms have no long 
range interaction. By coupling the ground and Rydberg 
states with a laser tuned A from the \g) — > \r) transition 
and with coupling strength il (the Rabi frequency) we 
mix the states \s) and |r). This technique [28] of dress- 
ing the atoms with the laser means that the trapped, 
ground state atoms acquire the characteristics of the Ry- 
dberg state, but in a controllable fashion. In particular, 
the coefficient C3 is replaced by an effective interaction 
coefficient which we write as C3 — > (il/2/S)C^. 

The atoms, must now be confined in a bilayer lattice. 
The "itinerant" layer represents e.g. the electrons in a 
condensed matter problem, and the "phonon" layer acts 
as a means of simulating the long range interactions be- 
tween electrons. For the purposes of this model, the itin- 
erant layer may be partially filled, but it is essential that 
tunneling between adjacent sites is allowed. It is equally 
important that the phonon layer is completely filled (with 
1 atom per site) and that the tunneling is forbidden - the 
Mott insulator phase [T7j . This is achieved by making the 
potential barrier in the itinerant layer smaller than in the 
phonon layer. A further complication arises because the 
atoms in the phonon layer must be set with low oscillat- 
ing frequencies, while the itinerant layer must be set up 



with high phonon frequencies. This can be achieved if 
the optical lattice in the phonon layer has a special con- 
figuration as seen in Fig. [1] This form may be achieved 
using painted potentials [29 or by superimposing higher 
order harmonics on a standard optical lattice. 

In the itinerant layer, fermions hop with amplitude 
t, and experience local Coulomb repulsion, U, as in a 
typical optical lattice experiment to simulate a Hubbard 
model. Lattice vibrations can be introduced to this sys- 
tem by displacing the atoms in the phonon layer. Atomic 
displacements do not affect the optical lattice, so the vi- 
brations of the atoms are momentum independent Ein- 
stein phonons with Hamiltonian, 
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and with polarization vectors £ fc „ = in orthogonal 
directions. Here uik.v = w v is the angular frequency of 
a phonon from mode v with momentum k and <P and 
d respectively represent creation and annihilation oper- 
ators for phonons. One interesting possibility in 2D is 
to distort the harmonic potential of the optical lattice to 
lift the degeneracy of the phonon modes. Also, the use 
of very distorted potentials allows the probing of anhar- 
monic effects. 

The introduction of small in-plane phonon displace- 
ments, tti, at site i causes the interaction between Ryd- 
berg states to become, Vki — — ft 2 /x 2 /4A 2 | — Ri 
«i| 3 , which can be Taylor expanded if the displacement 
is small, 
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The phonon fields can be quantized in the usual way by 
substituting u { = J2 ktl/ y/h/2N Muj k ^ k ^(d k ^e~ lk Ri + 

d kv e lk ' Ri ). Taking into account that all sites in the 
phonon layer are occupied, and utilizing the momentum 
independent modes, a multimode extended Hubbard- 
Holstein interaction is derived, 



H, 



cxt — Hol 



3^ 2 r> 2 

4A 2 



]T RirMd^+d^/Rjj (4) 



Here, Rij is a vector between an atom in the itinerant 
layer at site i, and an atom in the phonon layer at site 
j, rii is the number operator for fermions on site i and 
M is the mass of the atoms. The multimode nature is 
of particular interest because interaction with multiple 
phonon modes is difficult to simulate with current nu- 
merical techniques. 

A further simplification can be made by elongating the 
potentials in the phonon layer along the direction perpen- 
dicular to the planes, so that the Hamiltonian simplifies 
to the form, 
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FIG. 2. Comparison between Rydberg and Frohlich inter- 
actions for <I > (a;)/$(0) and various interplane distances. An 
excellent correspondence between the interactions can be seen 
for interplane distance 6 < a. Open points indicate the 
screened lattice Frohlich interaction and filled points inter- 
plane Rydberg interaction. Curves are shown in pairs cor- 
responding to the same $(a)/$(0), so that R sc — > oo cor- 
responds to b = 1.4218(0)o (circles), R sc = 1.4624(5) corre- 
sponds to b = a (triangles) and R ac — 0.33233(6) to b = a/2 
(squares) . 

with the interaction, 

9ij 4A 2 (6 2 + r 2 .) 5 / 2 ' l J 

where is the distance between the projection of sites 
i and j onto the same later and b is the interplane 
distance. For comparison, a standard Holstein model, 
where the local electron density couples to local opti- 
cal phonon modes [T3], has <fe = gSij. The Frohlich 
model describing the continuum limit of interactions be- 
tween electrons and strongly polarizable materials [3U] 
can be generalized to lattice models and has the form 
gij ~ exp(— rij/R sc )(b 2 + rf^) -3 / 2 where R sc is a screen- 
ing radius [3"T] (the latter are sometimes known as ex- 
tended Holstein interactions). 

The effective instantaneous interaction between 
fermions in the itinerant layer is <&(r) = J2m 9r,m,go,m- 
Fig. [2] shows a comparison between the shapes 
of Rydberg and lattice Frohlich effective interactions 
$(xai)/$(0) for various interplane distances and screen- 
ing radii. As would be done in the experiment, the near- 
neighbor interactions are matched by modifying b to get 
the closest possible correspondence to the Frohlich inter- 
action that is to be simulated. An excellent correspon- 
dence between the shapes of the interactions can be seen 
for interplane distance b < a. The inset shows the near- 
neighbor interaction for a range of b, for use as a reference 
to match other types of interaction. 

To demonstrate the scheme at relatively low phonon 
frequencies fiuj = t, we use quantum Monte Carlo to com- 
pute the properties of polarons and bipolarons in the bi- 
layer lattice and compare with results from the screened 
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FIG. 3. Polaron and bipolaron energies for different ranges 
of interaction within the Rydberg scheme proposed here, and 
for a screened Frohlich interaction where the near-neighbor 
interaction strengths have been matched. There is excellent 
agreement between the two schemes, demonstrating the effi- 
cacy of the simulator. Unless shown, error bars are smaller 
than the points. The lower plots show the residual between 
the Rydberg and Frohlich schemes. 



Frohlich interaction. Results from the simulations can be 
seen in Fig. [3j We also show the effects of the harmonic 
potentials used in conventional optical lattice schemes. 
For the low density system simulated, the harmonic po- 
tential makes little difference to the overall values of en- 
ergy and phonon occupation. A very close agreement is 
found between the quantum simulator and Frohlich sys- 
tem, further demonstrating the efficacy of our proposal. 
A key point here is that the quantum simulator is not lim- 
ited to the low density regime, and will be able to probe 
and tune key physics of unconventional superconductors 
at intermediate densities where there is most controversy 
about the mechanisms. 

Finally, we discuss some of the technical challenges 
which must be solved in order to realize this system. 
We expect that the simplest implementation would be 
in an optical lattice setup. So far, micron sized opti- 
cal traps have been demonstrated, together with direct 
imaging of single atoms on each lattice site [32] . Recently 
it has become possible to measure electronic dispersions 
[33] . Using optical lattices to create the bilayer system 
presents many challenges, in particular, the need for the 
two lattices to be as close a few wavelengths of the light 
A;, ideally 1.4A;. At the same time, the atoms in one of 
the lattices must be in the Mott insulating phase, whereas 



4 



the other must have appreciable tunneling. In present ex- 
perimental systems, tunneling between sites occurs when 
the lattice spacing A/ is around 0.5 /im, which is compa- 
rable to the diffraction limited spot size of any imaging 
system used to create the bilayer lattice. 

Alternatively, a bilayer optical lattice may be set up 
with painted potentials [291 , a powerful technique with 
a high level of control over the shapes of features in the 
optical lattice. Creating a bilayer painted lattice requires 
two focused horizontal sheets and beams focused between 
the sheets to form wider features in the phonon layer 
sheet than the itinerant layer sheet. While the painted 
potential scheme is harder to implement than traditional 
optical lattices, it avoids the need to deal with harmonic 
frequencies. 

In the bilayer system, the interesting physics are en- 
coded in the momentum distribution of the gas held in 
the itinerant layer. What we want to study is the cor- 
relation function C'(k, a; — k, — a). To extract this, the 
momentum distribution must be observable. In the ex- 
periments in Ref. 1341 this was accomplished using a time- 
of-flight method, where the momentum distribution of 
the roaming atoms is mapped onto the spatial location 
at the time of imaging. 

One further complication arises from employing a 
dressed atom-atom interaction. In the cases where the 
two atoms in the itinerant layer occupy a single lattice 
site, the dipole-dipole interaction may saturate to a con- 
stant value, independent of the atomic separation. Fol- 
lowing reference [25] we find that the separation R c at 
which this occurs is given as R c = {C^/2Ti\ Al) 1 / 3 (in 
this expression, C3 is the undressed C3 coefficient). This 
means working with relatively low principal quantum 
number n sa 20 where the undressed coefficient C3 ~ 10 
GHz flip?. In such a case, the atoms will interact as 
desired for lattices of period greater than 300 nm, but 
atoms sitting on top of each other will feel a spatially 
independent interaction. If necessary, the dimensionless 
interaction strength can be increased by lowering t by 
making deeper wells in the optical lattice. 

In this letter, we have shown how systems of cold Ry- 
dberg atoms in a bilayer can be used as a simulator for 
electron-phonon interactions in the presence of strong 
electronic correlation, of the type found in many un- 
conventional superconductors. We have carried out the 
mapping to an extended Hubbard-Holstein model and 
discussed the experimental requirements for constructing 
such a simulator. Using continuous-time quantum Monte 
Carlo simulations, we have shown that results from the 
simulator are comparable to a Hubbard-Frohlich model 
and determined further constraints on the mapping. Fur- 
thermore, we have discussed how the simulator can be 
achieved using contemporary techniques. The proposed 
quantum simulations system goes well beyond the pos- 
sibilities of previously proposed quantum simulators for 



the simulation of interactions with lattice vibrations. In 
particular, fermions can be simulated, arbitrary fermion 
density is achievable, and multiple phonon modes can 
easily be included. 
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